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We numerically investigate the structures of the near-plate temperature profiles close to the 
bottom and top plates of turbulent Rayleigh-Benard flow in a cylindrical sample at Rayleigh numbers 
Ra — 10** to Ra = 2 x 10^^ and Prandtl numbers Pr = 6.4 and Pr = 0.7 with the dynamical frame 
method [Q. Zhou and K.-Q. Xia, Rhys. Rev. Lett. 104, 104301 (2010)] thus extending previous 
results for quasi- 2-dimensional systems to 3D systems for the first time. The dynamical frame 
method shows that the measured temperature profiles in the spatially and temporally local frame 
are much closer to the temperature profile of a laminar, zero-pressure gradient boundary layer 
according to Pohlhausen than in the fixed reference frame. The deviation between the measured 
profiles in the dynamical reference frame and the laminar profiles increases with decreasing Pr, 
where the thermal BL is more exposed to the bulk fluctuations due to the thinner kinetic BL, and 
increasing Ra, where more plumes are passing the measurement location. 



One of the classical systems to study heat transport 
phenomena is turbulent Rayleigh-Benard (RB) convec- 
tion. In this system a layer of fluid is heated from below 
and cooled from above [H |2]. The dynamics and the 
global features of the system are strongly influenced by 
the properties of the flow in the boundary layers (BLs). 
Almost all theories that describe the heat transport in 
turbulent RB convection, from the early marginal stabil- 
ity theory [3] to the Shraiman & Siggia (SS) model [U [5] 
and to the Grossmann & Lohse (GL) theory [6], are es- 
sentially BL theories. Therefore, it is a key issue to fully 
understand BL profiles close to the horizontal plates. 

The key question in RB convection is how do the heat 
transport, expressed by the Nusselt number Nu, and the 
turbulence intensity, expressed by the Reynolds number 
Re, depend on the control parameters of the system? For 
given aspect ratio F = _D/L (D is the sample diameter 
and L its height) and given geometry, the control param- 
eters are the Rayleigh number Ra = PgAL^ / {kv) and 
the Prandtl number Pr = v/k. Here (3 is the thermal 
expansion coefficient, g the gravitational acceleration, A 
the temperature difference between the plates, and v and 
K, are the kinematic viscosity and thermal diffusivity, re- 
spectively. The GL theory has achieved good success in 
predicting Nu{Ra,Pr) and Re{Ra,Pr) [IIIIIIH]. 

Recently, the GL theory was successfully extended to 
the very large Ra number regime [3] (the so called ulti- 
mate range), in order to explain the experimentally ob- 
served multiple scaling regimes of the heat transfer [10] 
and to the rotating case to predict the heat transfer en- 
hancement due to rotation [11]. As the GL theory heav- 
ily builds on the assumption that the BL thickness of 
not yet turbulent BLs scales inversely proportional to 
the square root of the Re number according to Prandtl's 



1904 theory [T^], the degree of validity of Prandtl-Blasius 
BL [13] flow needs to be tested. Indeed, recent experi- 
ments [H] have shown that in non-rotating RB the BLs 
scaling behaves as in laminar flows. In addition, a com- 
parison of the mean bulk temperature calculated using 
the Prandtl-Blasius theory with that measured in both 
liquid and gaseous non-Oberbeck-Boussinesq RB convec- 
tion shows very good agreement |15H17j . Furthermore, 
the kinematic BL thickness evaluated by solving the lam- 
inar Prandtl-Blasius BL equations was found to agree 
well with that obtained in the direct numerical simulation 
(DNS) [18j . However, a comparison between experimen- 
tal velocity [T^] and numerical temperature [^ profiles 
obtained in the fixed reference frame (fixed with respect 
to the horizontal plates) with the respective classical pro- 
files shows significant deviations. 

Considerable progress on this issue was recently 
achieved by Zhou and Xia [21j who experimentally stud- 
ied the velocity BL for water (Pr = 4.3) with Particle Im- 
age Velocimetry (PIV). They argued that such deviations 
should be attributed to intermittent emissions of thermal 
plumes from the BLs and the corresponding temporal dy- 
namics of the BLs [51] [55] . This led to the study of the 
BL structures in dynamical reference frames, which fluc- 
tuate with the instantaneous BL thicknesses [21]. Zhou 
et al. [21] used PIV data to show that the velocity pro- 
files obtained in the dynamical reference frame are very 
close to the laminar Prandtl-Blasius profile and numer- 
ical data [221 m] to show that a good agreement with 
the temperature profile of a laminar, zero-pressure gra- 
dient boundary layer according to Pohlhausen [23j [24] , 
which we will call the "Pohlhausen" profile [5S] in the 
remainder of this paper is obtained . In addition, they 
showed that most of the time the rescaled instantaneous 



2 




FIG. 1. Comparison between temperature profiles obtained at the cell center axis near the bottom and top plate: dynamical 
9*{z*) (blue short dashed lines), fixed reference frame 6{z) (red solid lines), and the temperature profile of a laminar, zero- 
pressure gradient boundary layer according to Pohlhausen (black long dashed lines) for a (a) Pr — 6.4, Ra = 10*, and F = 1, 
(b) Pr = 6.4, Ra = 10*, and F = 1/2, (c) Pr = 0.7, Ra = 10*, and F = 1/2 (d) Pr = 0.7, Ra = 2 x 10^^ and F = 1/2. 



velocity and temperature profiles are also in agreement 
with the theoretical prediction. A dynamical BL rescal- 
ing method has thus been established, which extends the 
time- independent Prandtl-Blasius- Pohlhausen BL theory 
to the time-dependent case, in the sense that it holds lo- 
cally at every instant in the frame that fluctuates with 
the local instantaneous BL thickness. 

However, this dynamical rescaling method has up to 
now only been applied to experimental and numerical 
data obtained in (quasi-) 2D RB convection, where the 
orientation of the large scale circulation (LSC) is fixed. 
In addition, the experimental result is limited to Ra — 
3 X 10^^ . In this paper we will apply the dynamical rescal- 
ing method to previously obtained data from fully 3D 
DNS performed in a cylindrical sample. The questions 
we will address here are: Does the dynamical rescaling 
method also work for this 3D geometry, where more com- 
plex flow structures can form? Does the method work for 
the highest Ra number simulations, i.e. Ra = 2 x 10^^ 
we have available? And is a similar agreement obtained 
also for lower Pr? 

For the analysis presented here we use data from our 
previous simulations, i.e. the simulations of Stevens et al. 



[26l |27] for Pr = 0.7 in an aspect ratio F = 1/2 sample 
for Ra = 2x 10^ to Ra = 2 x 10^^, and the simulations of 
Stevens et al. HH] for Ra = 10^ Pr = 0.7 and Pr = 6.4, 
in F = 1 and F = 1/2 samples. All these simulations 
have been performed with a second order flnite differ- 
ence code, see Verzicco et al. [291 - I3T] for details. Dur- 
ing the simulations mentioned above we stored, besides 
a limited number of full 3D snapshots, a large number 
of vertical (and horizontal) snapshots over a period of 
30 to 150 dimensionless time units when time is scaled 
by L/y/ f3gAL. The vertical snapshots can be used to 
analyze the temperature profile in the BL at the center 
axis. This is a favorable location as the LSC is always 
passing the center axis, provided that the off-center (or 
sloshing) motion of the LSC is not too strong near the 
plates [321 [53]. In addition, restricting the analysis to 
the temperature field offers the benefit that the orienta- 
tion of the LSC just above the plates does not have to 
be determined, because the temperature profile should 
be independent of the azimuthal orientation of the LSC 
in the region around the center axis. This means that 
no difficult (and arbitrary) routines have to be used to 
determine the orientation and presence of the LSC. How- 
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FIG. 2. Comparison between temperature profiles obtained 
at the cell center near the bottom and top plate in the fixed 
reference frame (a) and the dynamical reference frame (b) for 
Pr = 0.7 in an aspect ratio F — 1/2 sample at Ra = 10** 
(red), Ra = 2 X 10^° (blue), Ra = 2 x 10^^ (dark green) 
and the temperature profile of a laminar, zero-pressure gra- 
dient boundary layer with Pr = 0.7 according to Pohlhausen 
(dashed line). Note that the difference between this theoret- 
ical profile and the numerical data increases with increasing 
Ra as indicated by the direction of the arrow. The insets 
in figure (a) and (b) compares the shape factor for the tem- 
perature profile of a laminar, zero-pressure gradient bound- 
ary layer with Pr = 0.7 according to Pohlhausen (horizontal 
dashed line) with the shape factor of the temperature profiles 
obtained in the fixed reference frame (panel a) and the dy- 
namical reference frame (panel b) . Note that the shape factor 
in the dynamical reference frame is closer to the theoretical 
prediction than in the fixed reference frame. 

ever, a drawback of this method is that we cannot deter- 
mine whether the LSC is passing the cell center at every 
time instance. As is shown by Zhou et al. [3T1 [53 [21] 
the agreement between the dynamically rescaled profile 
and the Prandtl-Blasius-Pohlhausen profile deteriorates 
when one considers a position where the LSC does not 
pass. Similar effects can influence our results, however 
preventing this is very difficult as it would require an ar- 
bitrary criterion to determine the presence of the LSC at 
a certain location. 



We determined the mean temperature profile in the 
dynamical reference frame by defining the instantaneous 
thermal BL thickness Sth{t) as the intersection between 
the temperature gradient at the plate and the tempera- 
ture of the first extremum value in the temperature pro- 
file [231 [21] ■ Even then we obtain that the mean bulk tem- 
perature obtained far away from the thermal BL is below 
9* . As discussed by Zhou et al. [23l [25 we attribute this 
to the emissions of thermal plumes. They propose to use 
the temperature at some position outside of the thermal 
BL, such as z/Xth = 3 or z^*^ = 3, as the asymptotic value 
for the BL, where Xth is the slope BL thickness obtained 
in the fixed reference frame and z'!^f^{t) = z/5th{t) is the 
rescaled vertical distance from the plate. In figure [l] we 
show a direct comparison between the various tempera- 
ture profiles, rescaled in the described way: the fixed ref- 
erence frame profiles 0(x, z)/9{x^ z — ZXth), the dynami- 
cal local profiles 0* [z^h) 1 0* {zlf^ — 3), where 9* indicates 
the temperature in the dynamical reference frame defined 
as 6l*(a;,z*J = {e{x,z = zl^5th{x,t),t)) [211 [21121] , and 
the Pohlhausen temperature profile. 

Figure [T^ and [i]d shows that for Ra = 10^ and Pr = 
6.4 the mean dynamically-rescaled profiles are very close 
to the Pohlhausen temperature profile profile for both 
F = 1 and F = 1/2 samples. In agreement with the 
results obtained by Zhou et al. [211 [231 [21] in (quasi) 
2D RB convection we find a significant preference of the 
dynamical frame based profiles: Around the thermal BL 
thickness the fixed reference frame based time averaged 
local profiles 6{z)/9{z = 3Xth) are all much lower than 
the Pohlhausen profile, while the dynamically rescaled 
instantaneous local profiles 9* {z'^fj / 9* {z'^i^ = 3) are much 
closer to the Pohlhausen profile. We emphasize that we 
observe this good agreement even though we do not filter 
out the moments in which the LSC is not passing the 
cell center. In addition, one has to realize that buoyancy 
driven instabilities generate a wall normal flow and break 
down the stratified structure of a laminar flow field in the 
BLs. Consequently, heat is no longer mainly conducted 
but also a significant fraction is convected throughout the 
BLs. This results in a deviation of the mean temperature 
profiles from the laminar BL predictions as seen in the 
figure. Another feature worthy of comment is that the 
results here show that the dynamic scaling method works 
also for F = 1/2 geometry, whereas all previous studies 
pT| [23l [24] only considered a F = 1 geometry. That the 
method still works in a F = 1/2 sample is remarkably 
given that the laminar BL theory has been developed for 
parallel flow over an inflnite fiat plate, whereas here in 
the aspect ratio F = 1/2 sample one can hardly find such 
regions of parallel flows at the top and bottom plates. 

In flgure [T]^ we show the results for Pr = 0.7 and 
Ra — 10* in a F = 1/2 sample. A comparison with the 
results for Pr = 6.4 at the same Ra reveals that for the 
lower Pr the dynamically rescaled profile differs more, 
but still less than in the fixed reference frame, from the 
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Pohlhausen temperature profile. The larger deviation at 
lower Pr can be explained by the fact that the viscous BL 
is thicker than the thermal BL for larger Pr and hence 
the thermal BL in this case is less influenced by the bulk 
turbulent flow than at low Pr. Figure [TJl shows that the 
difference between the dynamically rescaled profile and 
the Pohlhausen profile increases with increasing Ra. To 
quantify the deviations from the Pohlhausen profiles we 
determined the shape factor of the temperature profiles 
in the BLs close to the bottom and top plate, cf |Mj 



(1) 



where 5'^ and 5™ are the local displacement and momen- 
tum thicknesses of the profiles, respectively, defined as, 
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Here, Y = 0* oi 9 is the corresponding local temper- 
ature profile in the dynamical or fixed reference frame, 
and all ^-integrations are evaluated only over the range 
from z = to z = 3 instead of towards oo. Roughly 
speaking, the shape factor of a profile in general de- 
scribes how fast the profile approaches its asymptotic 
value, i.e. the larger the shape factor is, the faster the 
profile runs to its asymptotic level. The shape factor of 
the Pohlhausen temperature profile is H^^ « 2.59 for 
the present Pr = 0.7. Figure [2] shows that the shape 
factor is closer to the Pohlhausen value in the dynamical 
reference frame than in the fixed reference frame, which 
confirms that the profiles obtained in the dynamical ref- 
erence frame are indeed closer to the Pohlhausen profile. 
Here we emphasize that the peak in the PDF of the shape 
factor is even closer to the shape factor expected from the 
BL theory, e.g. for Ra = 2 x 10^° the peak in the PDF 
of the shape factor in the dynamical reference frame is 
aX H K, 2.50, whereas the average shape factor in the 
dynamical reference frame is « 2.42. Similar numbers 
are observed for the other Ra numbers. Furthermore, 
we see that the shape factor decreases for increasing Ra, 
which confirms that the difference between the dynam- 
ically rescaled profile and the Pohlhausen profile indeed 
increases with increasing Ra. We attribute this to the 
fact that the flow becomes more turbulent at lower Pr 
and higher Ra and a more turbulent flow leads to the 
formation of a larger number of plumes that pass the 
center axis. Because Zhou et al. [5T1 [231 Hi] have shown 
that passing plumes lead to a lower shape factor, it is not 
unexpected that the difference between the dynamically 
rescaled profile and the Pohlhausen prediction increases 
with Ra. In fact both for lower Pr and higher Ra the 
velocity field becomes more turbulent and advects the 



temperature differently than described by purely lami- 
nar theory which can leads to the larger deviations that 
are observed. We note, that even for the most turbulent 
case we considered here, i.e. Ra — 2x 10^^ and Pr = 0.7, 
we cannot identify any characteristic signature indicating 
that the BLs have become turbulent. In particular we do 
not observe logarithmic profiles, which would look signif- 
icantly different, i.e. very steep near the wall and much 
fiatter with extensions of order L away from the wall. 
Both these properties we do not observe in the simula- 
tions. This is in agreement with recent experimental [10] 
and theoretical [S] results, which show that the BLs be- 
come turbulent around 10^"^ < Ra < 5 x 10"'^''. However, 
we cannot exclude that the observed deviations in the BL 
profiles with respect to the laminar predictions are some 
preliminaries of what will come for significantly lower Pr 
or higher Ra. 

To summarize, we generalized the dynamical rescaling 
method of Zhou and Xia [211 [23l [24l for the first time 
to the data obtained from fully resolved 3D DNS of RB 
convection in a cylinder. In agreement with the results 
obtained in 2D we find that the method clearly reveals 
the temperature profile of a laminar, zero-pressure gradi- 
ent BL according to Pohlhausen in 3D flows. The method 
works best for relative low Ra and high Pr, where a rela- 
tively strong LSC is formed. For lower Pr and higher Ra 
the velocity field becomes more turbulent and advects 
the temperature differently than described by purely 
laminar theory. Therefore, the agreement between the 
Pohlhausen profile and the dynamically rescaled profile 
becomes less for lower Pr, which can be explained by the 
fact that the viscous BL is thicker than the thermal BL 
for high Pr and hence the thermal BL in this case is less 
influenced by the bulk turbulent flow than at low Pr. 
For high Ra the effect of passing plumes leads to a larger 
deviations with respect to the Pohlhausen temperature 
profile. 
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